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Abstract 

For non-relativistic Schrodinger equations the lowering of their degree by substitution 
\E'(r) — > F{r) = \E''(r)/\l/(r) is known to facilitate our understanding and use of their 
(incomplete, so called quasi-exact) solvability. We show that and how the radial 
Dirac relativistic equation may quasi-exactly be solved in similar spirit. 
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1 Introduction 



The two-component radial Dirac equation with energy E, mass M and centrifugal 
term f/(r) = fi;/r [jl| reads 



In a way similar to non-relativistic Schrodinger equation it proves exactly solvable for 
the electrostatic field of hydrogen atom V{r) = a/r accompanied, if necessary, by the 
auxiliary or external Lorentz scalar force W{r) = (3/r and by the possible central 



magnetic monopole charge Q m. n = + + 1 + 2\Q\) where i = 0, 1, . . .. 

Marginally, let us note that K = i+ l = Ois also admitted whenever > and 
g ^ (cf. ref. |] for details). 

Recently, Brihaye and Kosinski conjectured that a formal parallel between 
relativistic and non-relativistic quantum mechanics may be extended to many other 
models. Explicitly, they have demonstrated that the perturbation of hydrogen atom 
by the linear relativistic force W{r) ~ cur not only resembles its non-relativistic 
Coulomb plus linear plus quadratic analogue but also shares the incomplete, so called 
quasi-exact (QE) solvability with it. At certain exceptional energies and couplings, 
elementary bound states were obtained by non-numerical means in a way which 
complements the 25 years old non-relativistic result by A. Hautot 

In our present letter we intend to proceed one step further. Having in mind a deep 
non-relativistic connection between QE solvability and Riccati-Schrodinger equations 
1^, we shall formulate a parallel relativistic "order- lowering" idea and implement it 
in the technically slightly more difficult context of relativistic eq. (|l]). This will 
enable us to show that, in particular, the algebraic-equation approach of Brihaye 
and Kosinski just picks up a very specific portion of a much larger class of all the 
QE solvable Dirac equations. 
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2 New QE solutions: Explicit method 



Once we represent a wavefunction in non-relativistic quantum mechanics as an in- 
tegral Mr) = expg^^F{Od^ we get MrJ)/'^o{r,e) = [F{r)f + 9,F(r). This 
converts the radial differential Schrodinger equation to an equivalent first-order form 
V{r) = EQ — i{i + l)/r^+ [F{r)]'^ + drF{r). Such a Riccati-type re-arranged equation 
is nonlinear but may be re-interpreted as an explicit closed definition of a partially 
solvable potential in terms of any of its "tentative" wavefunctions. 

The latter point of view plays an important role in the understanding of non- 
relativistic QE systems P]. Mutatis mutandis, the relativistic pair of equations (|l]) 
may also define the QE solvable pairs of potentials W{r) and V{r). Indeed, assuming 
that the latter forces are responsible for the existence of any particular elementary 
wavefunction with components /o('") and go{r) at a particular mass Mq and energy 
Eq we may write 



W{r) 
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V(r) = -Eo + 
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On the basis of experience with non-relativistic models it is not too surprising that the 
mere correct threshold and asymptotic behaviour in tentative /o('") = pr'^exp(— Ar) 
and go{r) = qr^ exp(— A r) with normalization constants p and q already leads to the 
exactly and, incidentally, completely solvable model as mentioned above (cf. also 
1^). The next tentative elementary choice could mimick an unphysical singularity 
at negative r^, = — < 0, 

/o(r) = pr'^(l + /i r) exp(— A r), gQ{r) = qr^{l + hr) eyi]){—\r). (3) 

In terms of the same parameters e = ±1 and t G (— oo, oo) in the input ratio of 
norms p/q = s expt we get the same formula for energy E = —eX sinht and mass 
M = e X cosh t as above. Also both the Coulombic couplings remain the same, 

/ 

(4) 




- cosh t sinh t 
sinh t — cosh t 



e /i 

e K, 
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The only change emerges as a screening which enters the new and, by definition, QE 
solvable potentials 

a a 3 3 

V(r) = - + W(r) = -+ , , as = ehsmht, = -£ /i cosht. (5) 

r 1 + hr r 1 + hr 

Further states in such a model may be sought and studied, by the explicit algebraic 
method of ref. 0], in full analogy with semi-relativistic and non-relativistic QE 
solvable screened Coulomb potentials 0. 

3 QE solutions: Implicit method 

Our new QE model @ looks particularly simple after transition to the elementary 
integral representation of wavefunctions 

f{r) = e-^'™ , g{r) = e-'^'ir^i . (6) 

In non-relativistic setting, similar re-arrangement proved useful in computations 
as well as in the so called supersymmetric transformations of Hamiltonians . Also 
here, our postulate (|^) will lead to simplifications. Thus, with abbreviations Y{r) = 
[F{r) +G{r)]/2 and Z{r) = [F{r) — G{r)]/2, one of the components of Dirac equation 
lowers its order and becomes purely algebraic, 

[E + V{r)]^ + [r(r)]2 = [M + W{r)]^ + [f/(r) - Z{r)]\ (7) 

This enables us to parametrize, say, 

E + V{r) = R{r) cos A{r), M + W {r) = R{r) cos B (r) , (8) 

Y{r) = R{r) sin A{r), U{r) - Z{r) = R{r) sin B{r). (9) 

Such a transformation simplifies also the remaining Dirac equation 

M + W{r) - E -V{r) + [Y{r) + Z{r) -U{r)] exp2 F Z{^)d^ = (10) 

which, in terms of an abbreviation C{r) = —[A{r) + B{r)]/2, reads tanC(r) = 
exp 2 J^^^, Z{^)d^ or, in differential form, Z{r) = drC{r)/ sin2C(r). After an insertion 
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of such a definition of Z{r) in the second item of eq. we may ehminate the function 
R{r) and are left with a set of the closed and compact simultaneous definitions of 
both the wavefunctions and potentials in terms of an arbitrary initial choice of the 
auxiliary but practically unrestricted pair of functions A{r) and B{r). 



4 Re-parametrization 



Our above construction of a QE solvable system exhibits still a certain similarity 
to its non-relativistic Riccati-like version. Unfortunately, the parallel is incomplete. 
In particular, we cannot derive potentials immediately from the wavefunctions since 
both of them enter our formulae together with their derivatives. At the same time, 
our implicit QE-type solution of Dirac eq. (|^) still has to be made compatible with 
some overall physical requirements. Also its clearer physical interpretation is needed. 
For this purpose, let us make a further step. Recalling the standard Pauli matrices 

1 \ ( l\ 1 

^^y= \ > \ 

10/ 1-10/ I -1 



:iii 



let us re-write eq. (|ID as a manifestly real two-component problem 

{Idr + [M + W{r)]a^-t[E + V{r)]ay-U{r)a,} ipir) = (12) 

and pre-multiply it by the transposed two-component real spinors axi'ir) and iayipir) 
from the left. The resulting pair of relations 



(13) 



Z{r) — U{r) coshS(r) sinhS(r) E + V{r) 

—Y{r) sinhE;(r) coshS(r) M + W{r) 

with S(r) = 2 J^,^, Z{C,)dC, represents another, integral representation of our original 
Dirac bound-state problem. In it, the two-by-two matrix is easily invertible. 



E + V{r) 
M + W{r) 



coshS(r) — sinhS(r) 
— sinhS(r) coshS(r) 



Z{r) - U{r) 
-Y(r) 



(14) 



This induces the simplified ansatz 



E + V{r) = S{r)smhT{r), M + W{r) = S{r) coshT{r), 



(15) 
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Z{r) = U{r) + S{r) sinh[T(r) + S(r)], Y{r) = -S{r) cosh[r(r) + H(r)] (16) 

which, in particular, parametrizes all the QE solutions by the independent input 
functions S{r) and T{r). As already mentioned above, they must only be subject to 
the appropriate physical boundary conditions. We may conclude that the implicit 
relativistic implementation of the idea of QE solvability is as straightforward as its 
explicit nonrelativistic predecessor. 

5 Relativistic QE doublets 

After any change of our above "parametrization" point of view, technical complica- 
tions may re-emerge immediately. Even for the elementary and popular Coulomb -|- 
polynomial form of forces as suggested for further study of the QE solvability in ref. 
we immediately imagine that a seemingly trivial guarantee of their compatibility 
with our parameterizations leads in fact to a quite difficult algebraic problem. 

Another, mathematically easier and practically more important question is the 
possible existence and/or feasibility of constructions of the elementary QE multiplets. 
Indeed, in principle, after any choice of a QE wavefunction the resulting partially 
solvable potential may still remain compatible with another elementary bound state. 
Even in non-relativistic case such a physically useful requirement is mathematically 
non-trivial 0. There, it may still be characterized by the comparatively transparent 
condition 

dr[F,{r) - F^ir)] + [F,{r) - F2(r)][(Fi(r) + F^ir)] = - E,. (17) 

For the same QE potential (which was eliminated) this equation defines the superpo- 
sition Fs{r) = [Fi{r) + F2(r)]/2 in terms of the differences -Fd(r) = [Fi{r) — F2(r)]/2 
and 6 = E2 — El . This means that there exist very many nonrelativistic doublet 
partners Fi = Eg + Fii and F2 = Eg — E^ which are "numbered" by the choice of the 
virtually unconstrained functions Fd{r). 

In relativistic case, we have to proceed in similar vain. The elimination of the 
two energy-independent QE solvable potentials V{r) and W{r-) from the two versions 
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of eq. (T|) will describe the difference between the two right hand sides as an 
r— independent spinor proportional to 6. We get a relativistic differential-equation 
counterpart to eq. ([T7|) in terms of the auxiliary integrals a = a{r) = /^(^2 — Zi) 
and (3 = (3{r) = J^{Z2 + Zi) and their derivatives Z2{r-) = dra/2 + drl3/2 and 
Zi{r) = -dra/2 + dr (3/2, 

I 
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cosh(a + f3) 
sinh(a + /3) 



+ 



cosh 2a sinh 2a 
sinh 2a cosh 2a 




;i8) 



This formula may be read as an algebraic linear set of definitions of the two sums 
of exponents Yi{r) and Y2{r). The symmetry of our new equation with respect to 
the simultaneous double reflection j3{r) ^ —j3{r), U{r) ^ —U{r) and permutation 
Y'i(r) Y2{r) is one of the reasons why its closed solution is still unexpectedly 
compact, 



Y,(r) + Y2ir) 
YAr)-Y2(r) = 6 
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cosh j3{r) 
sinh a (r) 

sinh j3{r 



a\r) 



+ [p\r) - 2f/(r)] 



cosh a(r) 
sinh a (r) ' 

sinh afr) 



(19) 



(20) 



cosh Q;(r) cosh Q;(r) 

We may summarize that our "parameters" a{r) and /5(r) remain practically arbitrary 
functions. In the relativistic Dirac case there exist infinitely many QE bound-state 
doublets as well. Thus, the well known functional freedom of the implicit doublet 
solutions of the non-relativistic Schrodinger QE equations is shared by our present 
Dirac relativistic QE construction. 
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